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A time-variable, phase-dependent emission line in the X-ray 
spectrum of the isolated neutron star RX J0822-4300 
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ABSTRACT 

RX J0822^300 is the Central Compact Object associated with the Puppis A supernova rem- 
nant. Previous X-ray observations suggested RX J0822-4300 to be a young neutron star with a 
weak dipole field and a peculiar surface temperature distribution dominated by two antipodal 
spots with different temperatures and sizes. An emission line at 0.8 keV was also detected. We 
performed a very deep (130 ks) observation with XMM-Newton, which allowed us to study in 
detail the phase-resolved properties of RX J0822-4300. Our new data confirm the existence 
of a narrow spectral feature, best modelled as an emission line, only seen in the 'Soft' phase 
interval - when the cooler region is best aligned to the line of sight. Surprisingly, comparison 
of our recent observations to the older ones yields evidence for a variation in the emission 
line component, which can be modelled as a decrease in the central energy from ~0.80 keV 
in 2001 to ~0.73 keV in 2009-2010. The line could be generated via cyclotron scattering 
of thermal photons in an optically thin layer of gas, or - alternatively - it could originate in 
low-rate accretion by a debris disk. In any case, a variation in energy, pointing to a variation 
of the magnetic field in the line emitting region, cannot be easily accounted for. 

Key words: pulsars: general - stars: neutron - X-rays: individual: RX J0822-4300. 



1 INTRODUCTION 

Central Compact Objects (CCOs) in Supernova Remnants (SNRs) 
are a handful of about 10 point-like, thermally-emitting X-ray 
sources located close to the geometrical centres of non-plerionic 
supernova remnants, with no counterparts at any other wavelength. 
CC Os are supposed to be young, isolated, radio-quiet neutron stars 
(see lDe Luc 3120081 . for a review). 

While the first discovered CCO (the one in the R CW103 SNR) 
turned out to be a unique object jPe Luca et alj2006l) , results of X- 
ray timing on a sub-sample of sources with fast periodicities have 
recently set very useful constraints for a general picture of CCOs 
as a class. Analysis of multi-epoch XMM-Newton and Chandra 
observations of IE 1207.4-5209 inside 0296.5-^10.0 and CXOU 
J185238. 6-1-004020 at the centre of Kes 79 (P = 424 ms and 
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P = 105 ms, respective ly) yielded unambiguous evidenc e for very 
small period derivatives jHalpem & Gotthelfl2010l[201 ih . This im- 
plies, under standard rotating dipole assumptions, characteristic 
ages exceeding the age of the host supernova remnants by more 
than 3 orders of magnitude, as well as very small (<10" G) dipole 
magnetic fields. This points to an interpretation of such sources as 
'anti-magnetars' - weakly magnetised isolated neutron stars, born 
with a spin period almost identical to the currently observed one. 

Such a picture has been recently strengthened by the discovery 
of 1 12 ms pulsations from RX J0822-4300, the CCO in the Puppis 
A SNR (Gotthelf & Halpern 2009:), in two a rchival XMM-Newton 
datasets collected in 2001, previously used bv lHui&Becke3 ( l2006l) 
for a comprehensive spectral analysis of the neutron star and of 
the surrounding SNR. lOotthelf etallfeOlOl) , using a 2010 Chandra 
observation, set a 2cr upper limit of 3.5 x 10^^^ s s^^ to the period 
derivative, corresponding to a dipole magnetic field B < 2 x 10^^ 
G, and to a characteristic age Tc > 5 Myr, much larger than the age 
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of the host SNR (3.7 kvr Iwinkler etai]|l988l) . Thus, RXJ0822- 
4300 can be included in the anti-magnetar family. 

The X-ray emission propertie s of RXJ0822^300 are very 
peculiar. iGotthelf & HalpenJ hoOSh report on a 180° shift in the 
phase of the pulse peaks occurring abruptly at ~1.2 keV. This 
has been interpreted as due to the existence of two antipodal hot 
spots on the star surface, with two different temperatures ('warm' 
and 'hot') - indeed seen in the emission spectrum as two differ- 
ent blackbody components. The star rotation bringing into view or 
hiding such regions makes the observed spectrum to change from a 
soft phase (maximum alignment of the warm spot with the line of 
sight) to a hard phase (maximum alignmen t of the hot spot). Even 
more intriguing, iGotthelf & HalpemI l [2009l) report on the evidence 
for an emission line at ~0.8 keV, possibly associated with the warm 
spot. 

Here we report on a very deep XMM-Newton observation of 
RX J0822-4300, performed in 2009 and 2010, which allows us to 
study in detail the phase-resolved emission properties of the neu- 
tron star and to check for any time variability. 



2 OBSERVATIONS AND DATA REDUCTION 

Our study is based on a very deep (~130 ks) observation with 
XMM-Newton, originally scheduled to fit into a single satellite or- 
bit (revolution n. 1836), starting on 2009 December 17. However, 
launch of the Helios 2B spacecraft required support from XMM- 
Newton ground stations, which resulted in a ~50 ks data gap in the 
middle of the orbit. The observation was completed on 2010 April 
5. We also included in our study the two archival, shorter datasets, 
obtained on 2001 April 15 (~29 ks) a nd 2001 November 8 (~24 
ks), used bv lGotthelf & HalpemI dioO^ in their previous investiga- 
tions. A summary of the observations is reported in T able [T] 

W e focus on data obtained by the pn detector jStrtider et al.l 
l200lh of the European Photon Imaging Camera instrument (EPIC). 
The detector was operated in the small-window mode (time res- 
olution of 5.7 ms, field of view of 4'. 3 x 4'. 3). The thin optical 
filter was used in all observations. EPIC/pn Observation Data Files 
were processed with the most recent release of the XMM-Newton 
Science Analysis Software (SAS vl 1) using standard pipelines. 



3 TIMING ANALYSIS 

For our timing analysis, we selected the pn source events from a 
circular region of 30-arcsec radius, including only 1- and 2-pixel 
events (PATTERN to 4) with the default flag mask. Photon arrival 
times were converted to t he Solar system barycentre TBD time us- 
ing the Chandra position l lGotthelf & Halperrj|2009l). 

Before the analysis presented in IGotthelf & HalpemI l l2009l) . 
the pulsations in RX J0822-4300 eluded detection for many years, 
because of a phase shift of about half a cycle between the nearly 
sinusoidal profiles in the soft (< 1 .2 keV) and h ard bands. By com- 
puting Zl periodograms ( lBuccheriet"ai]|l983h . we found that the 
energy bands that maximise the pulsed signal are 0.46-1.17 keV 
and 1.25-5.12 keV. The resulting soft and hard pulse profiles of 
the individual XMM-Newton observations (Table[TJ were cross- 
correlated, shifted and summed to create two distinct pulse pro- 
file templates. Owing to the higher signal-to-noise of the 1.25-5.12 
keV profiles, we carried out the analysis in the hard band, but we 
checked that the results are fully consistent with those obtained 
in the soft band. The hard pulse profiles from temporal segments 
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Figure 1. Background-subtracted folded light curves for RX 10822^300 in 
the soft energy range (0.46-1.17 keV, upper panel) and in the hai'd energy 
range (1.25-5.12 keV. middle panel). The lower panel shows the hardness 
ratio (Hard/Soft), normalised to its average value. Phase intervals used for 
phase-resolved spectroscopy ai'e also marked. 



of the XMM-Newton observations were cross-coiTelated with the 
template to determine times of arrival at each epoc h. By means of 
a phase-fitting analysis (e.g. iDall'Osso et al.ll2003l') . we measured 
the periods given in Table [T](we also confirm the measurements by 
IGotthelf & Halpemll2009l for observations A and B). 

We attempted to obtain a full phase-connected timing solution 
(i.e. a solution that accounts for all the spin cycles of the pulsar) 
for the longest possible time. Unfortunately, it was not possible 
to univocally phase-connect the solution found for the two adja- 
cent datasets C and D to other observations (the uncertainty on the 
phase - propagated along the large time span to other observations 
- largely exceeds 1 cycle). 

A linear fit of the periods in Table[T] gives for the period 
derivative (9 ± 12) x 10^^^ s s~^, which translates into 3(t limits 
-2.7 X 10"^'' s s"^ < P < 4.5 X 10"^*^ s s~^ (in agreement wit h 
the limits recently published in IGotthelf. Perna & Halpenill2010l) . 
The period derivative at the Solar system barycentre results from 
the intrinsic pulsar spin-up/down plus the contributions due to 
any external g ravitational field and the pulsar proper motion (e.g. 
|Phinnevlll992h . I n the case of RX 1082 2^300, for an assumed dis- 
tance of 2.2 kpc ( iRevnoso et alj[T99^) . the Galactic contribution is 
negative and negligible (~— 2 x 10"'^" s s~^) and that produced by 
the proper motion (165 ± 25 mas vr~^ : Iwinkler & Petre|[2007l) is 
l-6lo.4 X 10"^'^ s s"\ While the total (Galactic -I- proper motion) 
contribution does not impact significantly on the current limits on 
the P of RX 10822^300, we note that it is larger than the period 
derivative measured for th e CCO in Kes 79 (~8.7 x 10"^* s s~^ 
iHalpem & Gotthelfll20I0l) . 

The phase of the pulse peak is energy dependent (see Fig- 
ure [T), with an offset of 0.44 ± 0.02 between the profile as seen 
at lower {E < 1.17 keV) and higher energy {E > 1.25 keV), the 
transition o ccurring quite abrup t ly at ~ 1.2 keV, consistent with the 
findings of IGotthelf & HalpenJ ( |2009|) . As shown in Figure |2| the 
pulsed fraction decreases from ~15% in the 0.4-0.6 keV energy 
range to <2% in the 1.1-1.3 keV range, then grows again to -^15% 
at_E> 2 keV, in good agreement with the model bv lGotthelf et al.l 
( l2010h . 
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Table 1. Journal of the XMM-Newton observations. " Mid-point of obser- 
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Figure 2. Energy dependence of the pulsed fraction (PF). The PF was eval- 
uated in overlapping 0.2 keV energy bins, incremented in 0.1 keV steps, as 
the ratio between the number of counts above the minimum and the total 
number of counts. Background has been subtracted. A clear trend is appar- 
ent, with a minimum in the 1.1-1.3 keV energy range. 



4 PHASE-RESOLVED SPECTROSCOPY 

Thorough phase-in tegrated spectroscopy of RXJ0822^300 has 
been published by iHui & Beckeil l l2006l) . We will focus here on 
phase-resolved spectroscopy. 

RX 10822-4300 lies in a very complex environment, which 
makes background subtraction a critical task. Using phase- 
integrated data, we evaluated an optimal selection of source events 
by maximising the signal to noise ratio in the 0.3-10 keV range as a 
function of the source extraction radius. The best choice turned out 
to be a circle of 17.5" centred at the Chandra position. We extracted 
the background spectrum from an annular region with radii of 28" 
and 35" centred on the source. Different source and/or background 
regions yield consistent results within ~2(t. 

We aligned the photon phases for each observation by cross- 
correlating the pulse profiles (see Sect. [3]l. We generated a com- 
bined event list fo r first-epoch observations (i.e. data collected in 
2001 and used by iGotthelf & HalDernll2009l) and one for second- 
epoch observations (i.e. our new data, collected in 2009-2010), in 
view of the large intercurring time span. We extracted four phase- 
resolved spectra from both first epoch and second epoch data, based 
on the intervals marked as 'Hard phase', 'Soft phase', 'Harden- 
ing' and 'Softening' in Fig. [T] Phase-resolved spectra were re- 
binned with at least 100 counts per bin and so that the energy 
resolution was not oversampled by more than a factor of 3. Re- 
sponse matrices and effective area files for each epoch were gener- 
ated by combining (weighted by exposure time) the corresponding 
files generated using the SAS tools RMFGEN and ARFGEN. Spec- 
tral modelling was performed with the XSPEC 12.6.0 package in 
the 0.3—4.0 keV energy range. To account for the effects of inter- 
stellar absorption, we used the TBABS model in XSPEC, with the 
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3. Residuals (in counts s~^) of two-blackbody fits to the Soft phase 
(see Fig. [T). A structure in the 0.6-0.9 keV range is seen in both 
although with different shape and intensity. 



abundances by Iwilms. Allen, & McCravl ( 120001) . We quote errors 
at 90% confidence level for a single interesting parameter unless 
otherwise specified. 

First, we repeated the exercise performed by 
iGotthelf &HalpenJ j2009i) , fitting a double blackbody model 
to our data. A simultaneous fit to the four phase-resolved spectra 
was performed for each epoch. The blackbody normalisations were 
allowed to vary as a function of the phase, while the temperatures 
and A^H were held fixed in order to constrain a single value for 
all phase ranges. This model yields a reduced-x^ of of 1.20 for 
298 degrees of freedom (d.o.f.) and of 1.23 for 389 d.o.f for the 
first and second epoch data, respectively. Although modulation of 
the emitting radii accounts for the bulk of the spectral variation, 
structured residuals in the 0.6-0.9 keV range are apparent in 
the Soft phase both in th e 2001 dataset (as already reported by 
IGotthelf & Halpemll2009h and in our deeper 2009-2010 dataset, 
which confirms the existence of a phase-dependent spectral feature. 

Very interestingly, while the phase-resolved best fit parame- 
ters do not change as a function of the epoch - they can be linked 
in a simultaneous two-epoch fit (more details below) - the devia- 
tion from the continuum in the soft phase has a somewhat different 
shape in 2001 with respect to 2009-2010 (see Fig.[3](, suggesting a 
possible time variability of the spectral feature. Indeed, such vari- 
ation is fully apparent when plotting together the two Soft phase 
spectra (see Fig. |4]l. To quantify the significance of the spectral 
change in a model-independent way, we compared the distributions 
of the source events' energies observed in the two epochs using the 
Kolmogorov-Smirnov test. The probability of a statistical fluctua- 
tion producing the apparent difference in the energy range where 
the feature is seen (~0.6-0.9 keV) turned out to be of 3 x 10~®. 

As a second step, in order to mod el the feature, we focused o n 
the two Soft phase spectra. Following IGotthelf & Halpenil ( l2009h . 
we added a Gaussian emission line to the two-blackbody model. 
Indeed, this yields a much better fit with no structured residuals in 
the 0.6-0.9 keV range. As expected, the line component varies as 
a function of the time, its central energy being higher in the first 
epoch (~0.80 keV) than in the second epoch (~0.73 keV). The 
significance of such line component was studied by calibrating the 
F-statistics using simula tions of the null model (the double black- 
body) as suggested by IProtassov et al. I l l2002h . Each epoch was 
studied individually. Running 10* Monte Carlo simulations, we es- 
timated that the significance of the line is greater than 99.99% in the 
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Figure 4. Upper panel: Soft phase spectra for the two epochs. First epoch: 
blue. Second epoch: green. The variation in the 0.6-0.9 keV range is appar- 
ent. The best-fit model, two blackbody components and a variable Gaussian 
emission line, is superimposed. The line component at both epochs is also 
shown. The continuum components do not change as a function of time. 
Lower panel: residuals to the best fit model. The model yields a very good 
description of the two spectra (xj = 1.02, 154 d.o.f.). 

second-epoch spectrum, and greater than 99.97% in the first epoch. 
Performing a simultaneous fit to the two spectra, the two-blackbody 
plus emission line model yields a reduced-x^ of 1.17 (155 d.o.f.) 
when all parameters (both continuum and line) are linked, while by 
allowing the line central energy to vary between first epoch and sec- 
ond epoch, the fit improves to a reduced-x^ of 1.02 (154 d.o.f.); see 
Fig- 111 The chance occurrence probability of such an improvement 
is ~3 X 10^^, as evaluated by an F-test. Such a result is confirmed 
by Monte Carlo simulations: assuming the best fit model of second 
epoch data, a line centroid as observed in first epoch could be ob- 
tained by chance with a probability as low as ~5 x 10~®. The fit 
does not improve significantly by allowing further line parameters 
to vary. There is no evidence for any variation in the continuum 
parameters between the two epochs. 

We also tried to model the two soft phase spectra using vari- 
able absorption features. Adding a single Gaussian absorption line 
(GABS in XSPEC) at a higher energy than the one required by the 
emission line model (~1.00 keV and ~0.89 keV in the first and 
second-epoch, respectively) yields a much worse fit than the emis- 
sion line model (reduced-x^ of 1.20 for 154 d.o.f.). Adding two 
lines at ~0.97 keV and at ~0.63 keV in the first epoch and at 
~0.88 keV and at ~0.55 keV in the second epoch still yields a 
remarkably worse description of the data (reduced ~ 1-15 for 
150 dof). As a further test, we tried lo w- 5 field neutron star atmo- 
spher e models for the continuum (NS A, IZavlin, Pavlov & Shibanovl 
1 19961 : NSATMOS, lHeinke et"al]|2006l) . As for the blackbody case, 
two components with different temperatures are needed. Addition 
of a variable Gaussian emission line is still favoured (x^ = 1-05 
for 154 dof) with respect to a single variable absorption line, as 
weU as to two variable absorption lines (x^ = 1-23 for 150 dof for 
the latter model). 

Then, the whole analysis was repeated for the other phase in- 
tervals (Hard, Softening and Hardening). No significant improve- 



ment in the fit was obtained by adding an emission line to the two- 
blackbody model (the same is true using absorption features). We 
assessed that, in each phase interval, the continuum did not change 
as a function of the epoch and that there are no systematic varia- 
tions between first and second-epoch in the 0.5-1 keV energy range 
(such results indicate that the long-term variability of the feature 
cannot have an instrumental origin). 

As a final step, we performed a simultaneous fit to all the 
spectra based on the results described above. We used the two- 
blackbody plus emission line model. The blackbody normalizations 
are phase-dependent, but not epoch-dependent; the line normalisa- 
tion is phase-dependent and its central energy is epoch-dependent. 
This yields a reduced x^ = 1-13 for 691 d.o.f. In such model, the 
A''h is (5.0 ±0.1) X 10^^ cm"^, the warm blackbody has a tem- 
perature fcTw = 265 ± 15 eV and a radius ranging from 2.27 km 
(Soft phase) to 2.04 km (Hard phase), the hot blackbody has a tem- 
perature kTn = 455 ± 20 eV and a radius ranging from 0.53 km 
(Soft phase) to 0.65 km (Hard phase). The line component is nar- 
row (o- < 40 eV). In 2001, the line energy is 0.80 ± 0.01 keV and 
the Equivalent Width (EW) ranges from ~53 eV in the Soft phase 
to < 10-15 eV in other intervals (where the line is not significant). 
In 2009-2010, the line energy is 0.73l[!:o2 keV and the EW ranges 
from ~45 eV (Soft phase) to < 12-18 eV in other intervals. The 
unabsorbed flux of the line between 0.3 and 10 keV is ~3% of the 
flux of the continuum in the same energy range. 



5 DISCUSSION 

Our multi-epoch XMM-Newton analysis shows a phase-dependent, 
time variable spectral feature, best modelled as an emission line 
with a variable central energy, in the X-ray spectrum of the 'anti- 
magnetar' candidate RX J0822^300. 

To put such a peculiar result in context, we first note that our 
observations confirm the picture of RX J0822-4300 as a weakly- 
magnetised neutron star. Indeed, we improved the upper limit on 
P, bringing the dipole component of the magnetic field down to 
< 2.3 X 10^^ G at 3(7 level. B ased on a larger statist ics, we also con- 
firm the geometric model bv lGotthelf et alj feOlOl) . explaining the 
phase-resolved thermal emission with two antipodal spots of dif- 
ferent temperature (compare e.g. our Fig.[2]to their Fig. 6). Lack of 
any measurable time variation in the continuum properties suggests 
that the warm and hot regions are intrinsic, persistent features in the 
thermal map of the star. To explain such a large surface anisotropy 
for RXJ0822-4300 (and for CCOs in general) we may consider 
that a large difference in intensity could exist between the internal 
and the external ma gnetic fields of the neutron star, as proposed by 
iTurolla et al.l ilOl ih to explain the properties of the low-magn etic 
field Soft Gamma Repeater SGR 0418-1-5729 teea et al.|[201ol) . In 
our case, an internal (toroidal -l- poloidal) field of a few 10^'^ G 
woul d be large enough to effectively charm el the heat flux from the 
core JGeppert, Kuker. & Pag3l2004l |200^ , but would not induce 
crustal fractures with consequent magnetar-like bursting activity. 

The most natural interpretation of the variable emission line 
is that of a cyclotron feature produced by electrons. If its cen- 
tral energy is associated to the fundamental cyclotron fre- 
quency, the magnetic field in the line emitting region would be 6- 
7 X 10^° (1 + z) G (where z ~ 0.25 is the gravitational redshift). 
This is quite compatible with the upper limit from the spin-down 
rate. The line is very narrow (a ^ 40 eV). If it is a cyclotron line 
then AE/E — AB/B and the relative variation of B over the 
emitting region (conservatively) needs to be ^ 10%. Thus, the line 
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should be produced in a very compact region. A variation of tlie 
central energy of the feature would require a change either in po- 
sition of the emitting plasma within a non- variable magnetic field, 
or in the intensity of the magnetic field itself. The lack of changes 
in the phase-resolved continuum emission rules out simpler, purely 
geometric explanations such as precession of the neutron star. 

To explain the generation of an emission line in the spec- 
trum of an isolated neutron star, the possibility of cyclotron scat- 
tering of surface thermal photons by a geometrically thin, optically 
thin layer of plasma could be considered. However, under simple 
assumptions (emission from the entire star surface; plane-parallel 
geometry; pure, conservative scattering), a scattering layer would 
produce an absorption line. One might invoke a spatially-limited 
scattering medium, possibly some distance away from the star sur- 
face. Photons coming from the part of the surface not covered by 
the layer could be scattered along the line-of-sight giving rise to 
an emission feature. The value of B derived from the line energy 
is somehow smaller than the upper limit on the surface field, so 
the line could indeed form at some height in the magnetosphere. A 
confined medium seems also to be required by the results of phase- 
resolved spectroscopy which shows that the emission line is seen 
mostly when the cooler spot is into view. Still, such a picture seems 
rather contrived, since the nature of the layer and the mechanism 
keeping the plasma suspended and confined in a compact blob re- 
mains to be understood. 

To ease the problem, an energy source unrelated to the 
surface thermal emission should be invoked to excite 
to higher Landau levels in t h e line emitting region. Indeed, 
iNelson. Salpeter. & WassermanI ( 1 19931) predicted that for neutron 
stars accreting at a low rate (Laccr < 10^* erg s^^), and endowed 
with magnetic fields of 10^^-10^"' G, accreting ions may lose en- 
ergy to atmospheric electrons via magnetic Coulomb collisions. 
Electrons, excited to high Landau levels, radiatively decay and part 
of the cyclotron photons are expected to escape producing an em is- 
sion line. According to lNelson. Salpeter. & WassermanI Tl993h . at 
B < 10^^ G, the fraction of the accretion-powered flux escaping in 
the line is expected to be very small (<c;5%), the largest part being 
reprocessed and emitted in a thermal continuum. Thus, one should 
postulate that the bulk of the X-ray luminosity of RX 10822^300 is 
accretion-powered, at variance with observations. It would require 
an accretion rate of ~ 2 x 10^'^ g s~^ implyin g - under standard 



relations for propeller spin-down iMenou et all 19991) - a P value 
m ore than 10 times larger than our upper l imit. Although the model 
bv lNelson. Salpeter. & WassermanI i ll 9931) does not fit to our case, 
low-level accretion of supernova fallback material (which cannot 
be ruled out, based on X-ray timing, as well as on the optical upper 
limits set by Mignani et al.2009) could play some role in generating 
an emission line in a low-B field atmosphere. A detailed investiga- 
tion of such possibility is beyond the scope of this paper. 

We will not go into further speculations about the line emit- 
ting mechanism. We stress that the evidence for time evolution of 
the spectral feature is model-independent and represents to date 
the first evidence for variability in an "anti-magnetar" candidate. 
Likely, such 'activity' is related to a variation of the magnetic field 
of the star. A ~10% decrease in ~8 yr seems vastly too steep to 
be attributed to the large-scale dipole field. Possibly, we are wit- 
nessing evolution of a localised multipole component, dominating 
close to the star surface. This would hint at the presence of a large 
internal field, as proposed to explain the anisotropic thermal map of 
the star. Precise X-ray timing, assessing the P and measuring the 
star dipole field, will add a crucial piece of information. Coupled to 
further sensitive phase-resolved spectroscopy to monitor spectral 



variability, this could help to solve the puzzles set by our results 
on RX 10822^300 which would have important implications for 
the understanding of the nature of CCOs and of their relations with 
other families of neutron stars. 
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